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ABSTRACT 

Galaxies moving through the intracluster medium (ICM) of a cluster of galaxies can lose gas via ram 
pressure stripping. This stripped gas forms a tail behind the galaxy which is potentially observable. 
In this paper, we carry out hydrodynamical simulations of a galaxy undergoing stripping with a focus 
on the gas properties in the wake and their observational signatures. We include radiative cooling 
in an adaptive hydrocode in order to investigate the impact of a clumpy, multi-phase interstellar 
medium. We find that including cooling results in a very different morphology for the gas in the tail, 
with a much wider range of temperatures and densities. The tail is significantly narrower in runs 
with radiative cooling, in agreement with observed wakes. In addition, we make detailed predictions 
of H I, Ha and X-ray emission for the wake, showing that we generally expect detectable H I and 
Ha signatures, but no observable X-ray emission (at least for our chosen ram-pressure strength and 
ICM conditions). We find that the relative strength of the Ha diagnostic depends somewhat on our 
adopted minimum temperature fioor (below which we set cooling to zero to mimic physics processes 
not included in the simulation). 

Subject headings: galaxies: clusters, galaxies: interactions, methods: N-body simulations 



1. INTRODUCTION 

Galaxies orbiting within a cluster may undergo a num- 
ber of interactions that are specific to dense environ- 
ments. These include interactions between intracluster 
gas and the galaxy, such as ram pressure stripping and 
starvation, interactions between pairs of galaxies, such 
as harassment, and interactions between the galaxy and 
the cluster gravitational potential, such as tidal strip- 
ping. The relative importance of these various mecha- 
nisms in transforming galaxy morphology, color, and gas 
content has been the subject of much debate, with large 
surveys helping to disentangle the various drivers (e.g. 
Moran et al 2007; van den Bosch et al. 2008). An alter- 
native way to gauge their relative importance, and to get 
a deeper understanding of the processes themselves, is to 
look for observational signs of the different interactions. 
This has been attempted recently in the context of ram 
pressure stripping (Vollmer & Hutchmeier 2007). One 
difficulty with this approach is relating the specific ob- 
servations to the underlying mechanism: how long does 
a given signature last and what does it tell us about the 
physical conditions in the galaxy and cluster? In this pa- 
per, we will use high-resolution numerical simulations of 
ram pressure stripping (RPS) to better understand the 
properties and fate of gas stripped during the encounter, 
and to predict and interpret observational signatures of 
RPS. 

The ram pressure of the ICM can strip gas from a 
galaxy moving at typical cluster velocities (Gunn & Gott 
1972), transforming a gas rich spiral to a gas poor SO. In 
addition, the galaxy ISM can be compressed, enhanc- 
ing the star formation rate (Byrd & Valtonen 1990; Fu- 
jita & Nagashima 1999). The ICM can also remove the 
loosely bound reservoir of gas that accretes onto a galaxy, 
thus slowly starving it of new fuel to form stars (Larson 
et al. 1980). In an earlier work (Tonnesen, Bryan & 
van Gorkom 2007), we studied cluster galaxies that had 



formed and evolved in a cosmological simulation, using 
an adaptive mesh refinement (AMR) code to highly re- 
solve a single cluster. The simulation included radiative 
cooling and star formation, and allowed us to predict the- 
oretically how infalling galaxies lose their gas. We found 
evidence for both starvation and stripping effects, show- 
ing that most galaxies that lost their gas did so without 
losing stellar mass, an indication that ram pressure strip- 
ping was an effective evolutionary mechanism. 

Observationally, a number of studies point to ram pres- 
sure stripping as a common environmental interaction. 
Spirals in the center of the Virgo cluster have smaller 
H I disks than stellar disks, indicating an interaction that 
does not affect the stellar component of galaxies (Cay- 
atte et al. 1990; Warmels 1988). Studies of H I deficiency 
have shown that galaxies in clusters have less neutral hy- 
drogen than their counterparts in the field (see the review 
by Haynes et al. 1984). Solanes et al. (2001) studied H I 
deficiency in a sample of 18 cluster regions, and found 
that H I deficiency decreases smoothly out to large pro- 
jected distances from cluster centers. However, the first 
sightings of head-tail structures in late- type galaxies were 
radio continuum observations of three galaxies in A1367 
(Gavazzi & Jaffe 1987). The discovery of these tails, with 
a maximum length of 30 kpc, was followed up in both H I 
and Ha (Sullivan et al. 1981; Gavazzi 1989; Dickey & 
Gavazzi 1991; Gavazzi et al. 1995). Observations of dis- 
placed H I within the optical disk in the direction of the 
radio continuum tails and H II regions at the edge of the 
"head" of the galaxy supported a ram pressure stripping 
scenario. 

In a recent, more sensitive and more highly resolved 
H I imaging study of Virgo, Chung et al. (2007) found 
a number of one-sided H I tails pointing away from 
the cluster center (tail lengths between 13 and 32 kpc). 
These galaxies are likely falling in for the first time and 
gas is already being removed at large projected distances 
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from the cluster center. In addition, detailed investiga- 
tions of a few individual galaxies using multiple wave- 
lengths have begun to unravel their probable histories 
(e.g. Growl et al. 2005; Chung et al. 2005). For exam- 
ple, observations of NGC 4522 indicate that the galaxy 
is undergoing ram pressure stripping, although it is out- 
side of the high density ICM (Kenney, van Gorkom, & 
Vollmer 2004). All of these observations indicate that 
ram pressure stripping does occur within the cluster en- 
vironment at a range of distances from the cluster center. 

Recent deep observations of clusters have revealed very 
long gas tails in H I (Oosterloo & van Gorkom 2005; 
Koopmann et al. 2008, Haynes et al. 2007). Oosterloo 
& van Gorkom find that the -120 kpc tail of NGG 4388 
in the Virgo cluster is well explained by ram pressure 
stripping. However, Koopmann et al. (2008) hesitate 
to attribute the tail of the Virgo pair NGG 4532/DDO 
137 to ram pressure stripping because it is more than 
1.5 times farther from the cluster center than any other 
H I tail in the Virgo cluster, and at 500 kpc is an order 
of magnitude longer than those found in the simulations 
of ram pressure stripping by Roediger & Briiggen (2008) 
and those observed by Ghung et al. (2007). Haynes et al. 
(2007) also observed a tail that is distant from the cluster 
center, and used the velocity distribution of the gas in the 
tail to conclude that the tail associated with NGG 4254 
was stripped due to galaxy harassment. The simulation 
by Due & Bournaud (2008) shows that a single fast flyby 
could have produced the tail. 

There have also been an increasing number of obser- 
vations of gas tails in Ra (Kenney et al. 2008; Gavazzi 
et al 2001; Yagi et al. 2007; Yoshida et al. 2004a,b; 
Yoshida et al 2008; Yoshida et al 2002; Sun, Donahue 
& Voit 2007). The Ha tails have filamentary structure 
and tend to be narrow. The suspected cause of the Ha 
emission varies from case to case. For example, Yoshida 
et al. (2004a) find that most of the emission near NGG 
4388 is from stripped gas that has been ionized by the 
galaxy's AGN. Gas heating leading to Ha emission could 
also be caused by thermal conduction from the IGM or 
turbulent shock heating as suggested by both Kenney et 
al (2008) and Yoshida et al. (2004a,b). In a different 
scenario, RB 199 is likely emitting in Ha due to star for- 
mation: it is probable that this galaxy has undergone a 
recent merger, making dense, star- forming gas easier to 
strip (Yoshida et al 2008). 

Tails have also been observed in X-rays, frequently op- 
posite a sharp edge in galactic X-ray emission or a bow 
shock (Sun et al. 2006; Kim et al. 2008; Irwin & Sarazin 
1996; Sun & Vikhlinin 2005; Machacek et al 2006; Wang 
et al. 2004). Sun & Vikhlinin (2005) conclude that 
the tail they observe is likely galactic gas that has been 
heated by supernovae explosions and then removed via 
Kelvin-Helmholtz stripping. Sun et al. (2006) interpret 
the tail of the small late-type galaxy ESQ 137-001 to be 
the cool stripped ISM mixed with the hot IGM. Wang 
et al. (2004) find a long, 88 kpc, tail from a disturbed 
disk-like galaxy. Like Sun et al. (2006), these authors 
conclude that it is a mixture of cool stripped gas from the 
galaxy and the ICM, largely because the X-ray emitting 
gas in the tail is cooler and denser than the surrounding 
IGM. Kim et al. (2008) find that the high metallicity of 
the tail from the large elliptical galaxy NGG 7619 indi- 
cates that it is gas from the galaxy, but as it is one of 



the dominant galaxies in its group, it is possible that the 
tail is actually sloshing of hot halo gas. 

Whether ram pressure stripping can cause these tails 
and how they survive in the IGM is unknown. It is un- 
clear how common long neutral tails are, as Vollmer & 
Hutchmeier (2007) examined the surrounding ^60 kpc 
around 5 H I deficient galaxies in Virgo, and did not find 
any excess stripped gas. Also, the relationship between 
the gas tails of different temperatures is poorly under- 
stood. In fact, three of the long gas tails mentioned above 
have been observed in multiple wavelengths. NGG 4388 
has a 120 kpc H I tail (Oosterloo & van Gorkom 2005), 
and an extended emission line region consisting of many 
faint gas filaments emitting in Ha close to the galactic 
disk (Yoshida et al. 2002; Yoshida et al. 2004a). As dis- 
cussed above, although this galaxy has an AGN that is 
heating stripped gas (Yoshida et al. 2004), gas heating 
leading to Ha emission could also be caused by ther- 
mal conduction from the IGM or turbulent shock heat- 
ing. Another galaxy with multiwavelength emission is 
ESQ 137-001 (Sun et al. 2006; Sun et al. 2007), with 
a tail detected in both Ha and X-ray. NGG 4438 also 
may have a very long tail observed in Ha (Kenney et 
al. 2008), while much smaller H I and CO tails have 
been observed (see Vollmer et. al. 2009 and references 
therein). Vollmer et al. (2009) found that some of the 
nearby Ha emitting gas is spatially coincident with GO 
gas, while some is not and is moving at higher velocities. 
There are three possible explanations for this difference 
in velocities: 1) differential stripping (less dense gas is 
stripped more efficiently from the galactic disk); 2) dif- 
ferential acceleration (less dense gas is accelerated more 
quickly by the stripping wind); or 3) a combination of 
differential stripping and differential acceleration. 

Simulations focusing on the tails of stripped galax- 
ies have found somewhat differing results. Roediger & 
Briiggen (2008), using an adaptive mesh refinement hy- 
drodynamics code, FLASH, find that for most ram pres- 
sures, observable tail lengths are about 40 kpc, with only 
three of their 31 tail length measurements having lengths 
over 80 kpc. However, these simulations did not include 
radiative cooling. Kapferer et al. (2009), using a smooth 
particle hydrodynamics code that includes radiative cool- 
ing and star formation, find that tails can survive with 
cool, dense clumps (T < 2 x 10^ K) out to well over 100 
kpc. 

In our previous work (Tonnesen & Bryan 2009, from 
now on TB09), we studied a detailed galaxy simulation 
that included radiative cooling in order to simulate a 
multiphase ISM. We found that radiative cooling led to 
a mixture of high and lower density gas in the disk, and 
the less dense gas was quickly stripped. The IGM wind 
was then able to stream through holes in the disk and 
ablate the surviving dense clouds. This is different from 
the scenario without cooling in which galactic gas is only 
stripped from the edges of the disk and any gas in the 
shadow of the galaxy is largely protected from the IGM 
wind. Although we focussed on the remaining gas disk in 
that paper, we noticed that the stripped tail was highly 
structured. We found that, although the densest gas 
clouds, with densities of molecular clouds, were not di- 
rectly stripped from the galaxy, there was still a lot of 
structure in the stripped tail, with relatively dense clouds 
surviving for hundreds of Myr. 
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In this paper, we run a set of high resolution simu- 
lations (about 40 pc resolution, which is small enough 
to marginally resolve giant molecular clouds) to under- 
stand how a multiphase ISM could affect the survival 
and structure of ram pressure stripped gas. It is im- 
portant to recognize that we do not attempt to include 
all of the physics involved in the ISM, focusing on how 
density fluctuations that are observed in the multiphase 
ISM of galaxies can affect gas tails. Following Roediger 
& Briiggen (2008), we focus on the stripped tails of gas, 
specifically examining the length, width, and substruc- 
ture of the tails (see also Roediger, Bruggen & Hoeft 
2006). Including radiative cooling allows us to estimate 
the density of the gas, and emission from H I, Ha, and X- 
ray gas separately. The questions that we highlight are: 
can cool gas tails survive in the ICM with turbulence and 
shock heating of the stripped gas? How can bright Ha 
and X-ray emission be induced in a tail of gas? Finally, 
to what velocities is stripped gas accelerated, and does 
differential stripping occur in gas disks? 

The paper is structured as follows. After a brief in- 
troduction to our methodology, we provide the general 
characteristics of our simulation (§2.1-3). We introduce 
the parameters of our specific simulations in §2.4. We 
then (§3) discuss our results, specifically focusing on the 
density, temperature, and velocity structure in the gas 
tail. In §4 we compare our results to observations. We 
discuss our choice of radiative cooling floor and resolu- 
tion in §5.1-2, and the impact of additional physics to 
our results in §5.3-5. Finally, we conclude in §6 with a 
summary of our results and predictions for observers. 

2. METHODOLOGY 

We use the adaptive mesh refinement (AMR) code 
Enzo. To follow the gas stripped from our simulated 
galaxy, we employ an adaptive mesh for solving the 
fluid equations including gravity (Bryan 1999; Norman & 
Bryan 1999; O'Shea et al. 2004). The code begins with 
a fixed set of static grids and automatically adds refined 
grids as required in order to resolve important features 
in the flow. 

Our simulated region is 311 kpc on a side with a root 
grid resolution of 128^ cells. We allow an additional 6 
levels of refinement for our runs that include radiative 
cooling, for a smallest cell size of 38 pc. We refine our 
simulation based on the local gas mass, such that a cell 
was flagged for refinement whenever it contained more 
than 4892 Mq (our standard, best mass resolution), or 
41745 Mq (our lower mass resolution, used in TB09). We 
found that both parameters refined most of the galactic 
disk to 38 pc resolution by the time the wind hit, al- 
though the disk fragmentation was better resolved in the 
higher resolution run, resulting in smaller mass clumps. 
The high mass resolution run also refined much of the 
wake to a spatial resolution of about 76 kpc, while the 
run with lower mass resolution generally resulted in a 
factor of 2-4 lower spatial resolution in the wake. 

The simulation includes radiative cooling using the 
Sarazin & White (1987) coohng curve extended to low 
temperatures as described in Tasker & Bryan (2006). To 
mimic effects that we do not model directly (such as tur- 
bulence on scales below the grid scale, UV heating, mag- 
netic field support, or cosmic rays), we cut off the cool- 
ing curve at a minimum temperature Tmin so that the 



TABLE 1 

Galaxy Stellar and Dark Matter Constants 



Variable 


Value 


M* 


1 X 10^1 M0 




4 kpc 


b* 


0.25 kpc 


^bulge 


1 X 10^° M© 


^bulge 


0.4 kpc 


^DM 


23 kpc 


PDM 


3.8 X 10-2^ g cm-3 



cooling rate is zero below this temperature. In the sim- 
ulations described here we use either Tmin = 8000 K, or 
^min = 300 K. Both allow gas to cool below the thresh- 
old for neutral H formation. In TB09, we found that 
the minimum temperature affected the range of masses 
and sizes of clouds forming in the disk. This, in turn, 
resulted in changes in the timescale for gas loss from the 
disk, although the total amount of gas lost was similar. 
In this paper, as we will show in more detail later, we find 
that the structure of the wake also depends somewhat on 



2.1. The Galaxy 

Our galaxy is placed at a position corresponding to 
(155.5,155.5,68.42) kpc from the corner of our cubical 
311 kpc computational volume, so that we can follow the 
stripped gas for more than 200 kpc. The galaxy remains 
stationary throughout the runs. The ICM wind flows 
along the z-axis in the positive direction, with the lower 
z boundary set for inflow and upper z boundary set as 
outflow. The X and y boundaries are set to outflow in all 
three of our high resolution cases (8000K, 300K, and no 
cooling) and reflecting in the lower resolution cases. The 
different boundary conditions in the two runs had only a 
very minimal impact, since the box is so wide (311 kpc 
wide). 

We chose to model a massive spiral galaxy with a flat 
rotation curve of 200 km s~^. It consists of a gas disk 
that is followed using the adaptive mesh refinement al- 
gorithm (including self-gravity of the gas), as well as the 
static potentials of the stellar disk, stellar bulge, and dark 
matter halo. We directly follow Roediger & Bruggen 
(2006) in our modeling of the stellar and dark matter 
potential and gas disk. In particular, we model the stel- 
lar disk using a Plummer-Kuzmin disk (see Miyamoto & 
Nagai 1975), the stellar bulge using a spherical Hernquist 
profile (Hernquist 1993), and the dark matter halo using 
the spherical model of Burkert (1995). This dark matter 
halo model is compatible with observed rotation curves 
(Burkert 1995; Trachternach et al. 2008). The equation 
for the analytic potential is in Mori & Burkert (2000). 

The gas is described as a softened exponential disk: 

p{R, z) = ^ 0.5^sech (^) sech (j^) (1) 

Given this gas density distribution in the disk, the initial 
gas temperature and pressure are calculated to maintain 
the disk in hydrostatic equilibrium with the surrounding 
ICM in the z direction. The gas disk's rotational velocity 
is set so that the combination of the centrifugal force and 
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TABLE 2 
Gas Disk Constants 



MgaS 


1 X 10^0 Mo 


^gas 


7 kpc 




0.4 kpc 



the pressure gradient of the disk balances the radial grav- 
itational force. We taper the gas disk smoothly by multi- 
plying the gas density distribution by 0.5(1 +cos(7r(R— 20 
kpc)/26 kpc)) for 20 kpc < R < 26 kpc. See our galaxy 
parameters in Tables [T] and [21 

To identify gas that has been stripped from the galaxy 
we also follow a passive tracer which is initially set to 
1.0 inside the galaxy and 10~^^ outside. In the following 
analysis, we will use a minimum tracer fraction of 25% 
to find gas stripped from the galaxy (our conclusions do 
not change when we use 10%). 

2.2. ICM Conditions 

The galaxy initially evolves in a static, high-pressure 
ICM with p = 9.152 x 10"^^ g cm-^ and T = 4.15 x 
10^ K, to allow cool, dense gas to form in the galaxy. 
This naturally generates a multiphase ISM (see Tasker 
& Bryan 2006 and TB09 for more discussion of the ISM 
properties). 

After 155 Myrs, we reset the boundary conditions to 
generate a constant ICM infiow along the inner z-axis, 
which is always face-on to the galaxy. We choose to 
study the highest of the three ram pressure strengths ex- 
amined in TB09 (this case results in the largest mass in 
the wake). We choose Pram = P'^^cm = 6.4x10"^^ dynes 
cm~^. For the corresponding ICM physical conditions, 
p, T, and '^icm, we use the results from our earlier cluster 
simulation to find the mean density, velocity, and tem- 
perature encountered by infalling galaxies for a typical 
cluster mass, given the ram pressures value (Tonnesen, 
Bryan & van Gorkom 2007). See TB09 for other details 
regarding the numerical setup (we use case PHRCW in 
TB09). Because we are interested primarily in the im- 
pact of radiative cooling on a gas stripped tail, we only 
consider a case with a face-on wind, again because this 
will result in the most gas in the tail. 

2.3. Suite of Simulations 

We run two main simulations of a ram-pressure 
stripped galaxy including radiative cooling using a 
temperature-dependant cooling curve from Sarazin & 
White (1987) extended to either T^in = 8000 K or 
Pmin = 300 K, as described earlier. In addition, we run 
a simulation without radiative cooling. In this way, we 
can investigate the impact of cooling to various temper- 
ature fioors. These three runs all have the highest mass 
resolution and provide the main results of this paper. 

All simulations have the same maximum refinement 
level of / = 6 (i.e. the same best spatial resolution of 38 
kpc), but we can compare both of our radiatively cooled 
runs to ones with a factor of 8 lower mass resolution (as 
described earlier). This allows us to investigate the effect 
of numerical resolution. 

A final comparison simulation has the same refinement 
criteria as our highest mass resolution runs, but calcu- 



lates cooling rates by following the non-equilibrium ion- 
ization fractions of hydrogen and helium, and directly 
computing their cooling rates. We do not include metals 
in this cooling calculation, which leads to slower cooling 
of hot ISM gas. This in turn results in a stripping rate 
which lies between our radiative cooling runs and a run 
without radiative cooling (largely because of the smaller 
amount of fragmentation in the disk). This comparison 
run was primarily done to allow us to check our simple 
calculation of Ha emission - described in the next section 
- by using the hydrogen ionization state as computed by 
the code. 

2.4. Projections 

Enzo outputs the density and temperature of the gas in 
each cell. To transform these values into H I column den- 
sity and Ra intensity, we used Cloudy, version 08.00 of 
the code last described by Ferland et al. (1998). Using a 
grid of temperatures and densities, we calculated the hy- 
drogen neutral fraction and Ra emissivity. In the Cloudy 
calculation, we included CMB radiation, the cosmic ray 
background, bremsstrahlung radiation from the ICM and 
the 2005 version of the Haardt & Madau (2001) z = 
metagalactic continuum, as implemented by Cloudy. We 
found that including the local interstellar radiation field 
emission resulted in lower amounts of neutral gas. Since 
much of our gas is very distant from the galaxy, we de- 
cided not to include this radiation. We also found that 
removing bremsstrahlung radiation did not change any 
of the values we considered. 

We chose to calculate the neutral fraction and Ha emis- 
sivity for a thin plane-parallel gas cloud of width 100 pc. 
We selected this width because it loosely corresponds to 
the cell size of most of the gas in the highly resolved tails, 
and accounts approximately for radiative transfer effects. 
If we assumed the radiative thin limit, it would decrease 
the amount of H I we predict, and increase the Ha emis- 
sion for dense, low-temperature gas. Ideally, we would 
include the radiation field with radiative transfer directly 
in the simulation, but this is not yet feasible (and only 
has a mild impact on the dynamics); instead we post- 
process these results to get reasonable predictions for the 
ionization fraction and Ha emissivity (see Furlanetto et 
al. 2005 for a discussion of various approaches in the 
context of Lya emission). We do compare our projec- 
tions with ones made assuming clouds sizes of 10 pc and 
1 kpc. 

To create X-ray surface brightness projections, we use 
a spectral lookup table for low density gas with a metal- 
licity of 0.3 solar computed using a Raymond- Smith code 
(Raymond & Smith 1977), as updated in XSPEC (Ar- 
naud 1996). The X-ray band we use is 0.5 keV to 2.0 
keV, following Sun et al. (2006). 

3. TAIL ATTRIBUTES 

In this section, we examine the physical characteristics 
of the wakes, focusing on the wake morphology, phase di- 
agrams and the velocity of the stripped gas. Then, in the 
following section, we turn to observational diagnostics. 

3.1. Morphology 

We first consider the morphological characteristics of 
our simulations. See Figure [T] to compare the surface 
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Fig. 1. — Surface density of all gas 500 Myr after the wind has hit the galaxy for our runs with radiative cooling to Tmin = 300 K (left), 
^min = 8000 K (center) and for the comparison run with no radiative cooling (right). The images are shown with a logarithmic stretch. 
Note the large amount of structure, the length, and the lack of flaring in the tail with radiative cooling, compared to the tail without 
radiative cooling. The radiatively cooled tails are cut off in this projection 212 kpc from the galaxy. 



density of the T^in = 300 K (left), T^in = 8000 K (cen- 
ter) simulations and the run without radiative cooling 
(right) at a time 500 Myr after the wind has hit the 
galaxies. Clearly there are important differences in the 
morphology of the resulting tails, although, remarkably, 
the amount of mass lost by this time differs by only 20% 
between Tmin = 8000 K and the run without radiative 
cooling. 

We find there are three primary changes that radiative 
cooling makes in the tail morphology. First, because of 
the formation of clouds in the disk, the runs with cooling 
show more structure in the tail. As we will see below, 
there is a wide range of densities and temperatures, with 
lots of small pressure-confined clouds. This is entirely 
missing in the no-cooling run and makes a crucial dif- 
ference in some of the observational diagnostics. This 
result is largely independent of Tmin, as long as cooling 
can generate a clumpy ISM. 

Second, the tails in the cooling runs are more extended 
at a given time, reaching more than 200 kpc past the 
galaxy within 500 Myr of the wind hitting the galaxy. 
In the run without radiative cooling, some gas has gone 
as far as 150 kpc, but the majority of the tail is closer 
than 125 kpc from the galactic disk. In all of our simula- 



tions, stripped gas begins moving away from the galaxy 
as a group, however, in the cases with radiative cooling, 
with time the material is stretched out as the velocities 
vary (we address the velocity distribution in more detail 
later). 

Finally, the stripped gas in the radiative cooling runs 
expands only slightly in the transverse direction as it 
moves with the ICM wind: in the Tmin = 8000 K run, 200 
kpc from the disk the tail width is about 80 kpc (500 Myr 
after the wind has hit the galaxy). The Tmin = 300 K 
run is even narrower, with a tail width of 40 kpc 200 kpc 
downstream of the disk. This is much less flaring than 
seen in in our no-cooling run: 500 Myr after the wind 
has hit the galaxy, the tail is 85 kpc wide at distance 125 
kpc downstream from the disk. Although this is total gas 
surface density, which cannot be observed, it is notable 
that most observed tails are long and narrow, and do not 
seem to flare (for example, see the images in Chung et 
al. 2007). 

The no-cooling results are in qualitative agreement 
with the simulations of Roediger & Briiggen (2008), who 
also find the tails to have a decreasing density profile 
with distance from the galaxy. We find that including 
radiative cooling permits higher density and lower gas 
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1 0- [ , , , 500,..Myr j ^ ^ ^ ^ ^ ] ^ ^ ^ ^ ^ ; 

10^ 10^ lO'* 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 10^ 
T (K) T (K) T (K) 

Fig. 2. — Contour plots showing the mass in gas at different densities and temperatures for our Tmin = 300 K run on the left, Tmin = 8000 
K in the middle, and the comparison run without radiative cooling on the right. The contours are log mass. For a gas cell to be included 
in the contour plot, it must have at least 25% of its mass originating from the galaxy. The line denotes the surrounding ICM pressure. 
Note that at lower temperatures the gas is below the ICM pressure due to rapid cooling (see text). Note also that the range of pressures 
in the simulations with cooling are much larger than the comparison run without. See section [3. 2 1 for discussion. 



to form in the disk and hence to be stripped, leading to 
high density gas throughout the tail. Our surface den- 
sity projection grossly resembles that of Kapferer et al. 
(2009), who also include radiative cooling in their sim- 
ulations, in that we both have elongated gas clumps in 
the wind out to large distances above the galactic disk. 
However, at less than a kiloparsec in diameter, most of 
our clumps are much smaller (about 10% in linear size) 
and much less massive than the clumps in the Kapferer 
et al. (2009) simulation. 

3.2. Gas Temperature and Density 

In Figure [21 we show the mass- weighted distribution of 
density and temperature for gas in the wake that origi- 
nated from the galaxy. The plots include all of the gas 
located between 10 kpc and 240 kpc above the disk that 
has at least 25% of the gas originating in the galaxy. We 



show these plots for three different times after the wind 
has hit the galaxy: 100 Myr, 250 Myr, and 500 Myr for 
the runs with radiative cooling to 300 K (left), with ra- 
diative cooling to 8000 K (middle) and without cooling 
(right). The solid line shows the background ICM pres- 
sure of 1.76 X 10~^^ dynes cm~^. 

Once again, we see a clear difference between the runs 
with and without radiative cooling. To begin, we concen- 
trate on a comparison between the no-cooling run and 
the simulation with cooling to Tmin = 8000 K. We ex- 
plore three distinct changes in the phase diagram. 

First, there is the simple fact that gas exists at much 
higher density and lower temperature in the runs with 
cooling, even at constant pressure. This is mostly due 
to the formation of dense clouds in the galaxy ISM be- 
fore stripping (see also TB09), but it is clear that cooling 
plays a role even at late times: in the no-cooling run, the 
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Fig. 3. — Zoom-in of a region of the wake from the run with 
coohng to 8000 K near the galactic plane 180 Myr after the wind 
has hit the disk. The greyscale shows temperature (lowest tem- 
perature is black), while the red contours indicate pressure. The 
outermost contours have the largest pressure. The slice is 6 kpc in 
depth, and 20 kpc x 37 kpc in size. Note that the lowest pressure 
contours surround the coolest gas clouds (in black). See section [321 
for a discussion. 



distribution becomes narrower as time goes on, as the gas 
is heated and mixed, while, if anything, the coohng dis- 
tributions become more extended. This higher- density 
gas shows up clearly in the tail morphology as shown in 
Figure [H 

Second, the cooling runs show a wider spread of density 
and pressure at fixed temperature; this is true even at 
large temperatures, but appears to be most noticeable 
at lower temperatures. Turbulence can cause a spread 
in pressure throughout the gas in the wind, as can be 
seen both in the run without radiative cooling, and at 
higher temperatures in our radiatively cooled runs, where 
cooling is inefficient. 

The final distinctive difference is a clear break in the 
distribution at T ^ 10^ K, with gas below this temper- 
ature falling off the constant pressure line. This corre- 
sponds to the peak of the cooling curve, when the cooling 
time becomes very short. In cooler clumps of stripped gas 
there is an interplay between cooling, which lowers the 
pressure of the gas, and adiabatic compression or shock 
heating from the surrounding ICM. At T ~ 10^ K, radia- 
tive cooling decreases the temperature more quickly than 
compression can increase the pressure, resulting in most 
of the cool gas lying below the pressure of the ICM. The 
peak in the gas pressure in the Tmin = 8000 K run oc- 
curs where radiative cooling is turned off, which suddenly 
allows the gas to be recompressed to the ICM pressure. 

To investigate this process in more detail, we look 
closely at a slice from our Tmin = 8000 K simulation 
in Figure [3l in which grayscale contours of temperature 
(with black being the lowest temperatures) are plotted 
with red contours of pressure, with the outermost con- 
tours denoting the highest pressure. This 6 kpc thick 
slice is 20 kpc x 37 kpc, and is taken 180 Myr after the 
wind has hit the galaxy. As expected from the above ar- 
gument, the lowest temperature regions are also the lower 
pressure regions, indicating that cooling lowers pressure 
faster than adiabatic compression can increase it. 

A simple timescale argument shows that the cooling 



and compression times are similar. For example, at T = 
10^ K, the cooling time for a clump density of 0.1 cm~^ 
is about the same as the sound crossing time across the 
clump, given a size of 1 kpc. As we find that many of 
our clouds are smaller than this size, compression must 
be acting more slowly than estimated using this simple 
calculation. 
The 

^min — 300 K run is similar but shows that gas 
can continue to cool below 8000 K in dense regions. The 
return to high pressures can be seen at this lower Tmin, 
otherwise the gas is in a rough pressure equilibrium below 
T ~ 10^ K, but at a pressure approximately an order of 
magnitude below the ICM pressure. 

3.3. Velocity Distribution 

In this section we examine the velocity structure of the 
stripped gas in our three main runs, from left to right: 
Tmin = 300 K, Tmin = 8000 K, and the no-cooling run. 
Figure [H shows the gas z- velocity, parallel to the ICM 
wind at times 100 Myr, 250 Myr, and 500 Myr after the 
wind has hit the galaxy. Again we are only considering 
gas that originated from the disk (a 25% minimum tracer 
fraction). We plot contours of gas mass as a function of 
velocity and distance above the disk. In our runs with 
radiative cooling, we find that the majority of stripped 
gas is never accelerated to the ICM wind speed of 1413 
km/s. In fact, 500 Myr after the wind has hit the galaxy, 
the main contour of gas asymptotes at about 900 km/s 
from about 150 kpc above the disk in both cooling runs. 

Our comparison run without radiative cooling evolves 
very differently. All three cases look similar in velocity 
space 100 Myr after the wind has hit the galaxy, but 
by 250 Myr - in the no cooling case - the main mass 
contour of stripped gas is accelerated to very close to the 
ICM speed (1413 km/s) at a height 100 kpc above the 
disk, with a narrow velocity width of only a ~ 250 km/s. 
The radiative cooling cases accelerate more slowly with 
distance above the disk, with none of the main contour 
of stripped gas at the ICM velocity at 250 Myr. By 
500 Myr after the wind has hit the disk, the radiative 
cooling cases show a narrow range in their tail velocity 
distributions, but the case without cooling has a very 
wide range of velocities, with stripped gas moving at both 
the ICM velocity and falling back towards the disk (at 
a height of 150 kpc above the galaxy). The width of 
the wake in the velocity distribution is a measure of the 
strength of turbulence seen in the wake. At late times, 
the radiatively cooled cases have a characteristic width 
of cr ~ 700 km/s, while the no-cooling case seems to be 
more turbulent, with a characteristic width of a ~ 1000 
km/s. 

To understand this difference, we turn back to the dif- 
ferent morphology of the tails. The important difference 
between our cooling and no-cooling runs is that radia- 
tive cooling forms much more structure in the tail. These 
dense clumps may be more difficult to accelerate and less 
affected by turbulence, so will not shift laterally into the 
shadow of the disk. In the no-cooling run, the large vari- 
ation in the wind velocity has to do with the large eddies 
that form from the interaction between the stripped gas 
and ICM, as seen in Roediger & Briiggen (2007) (their 
Figure 6). These eddies move the gas into the shadow 
of the disk, where the gas can fall back onto the galaxy. 
Also, both this work and Roediger & Briiggen (2008) find 
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Fig. 4. — Contour plots of gas mass at different velocities parallel to the wind (perpendicular to the disk) for the Tmin = 300 K, 
^min = 8000 K, and no cooling cases (going left to right). The contours are plotted as a function of height above the disk. Again, only gas 
cells with at least 25% of their m ass o riginating from the galaxy are counted. Very little of the disk gas is accelerated to the ICM wind 
velocity of 1413 km/s. See section [231 for a discussion. 



that some stripped gas is accelerated to the ICM veloc- 
ity about 100 kpc downwind of the galaxy (this is true 
even though Roediger & Bruggen (2008) vary their wind 
velocity while ours remains constant at 1413 km/s). 

In addition to looking at the velocity structure of all 
stripped gas, we can see how different density gas is af- 
fected by the ICM wind, looking for differential stripping 
and acceleration. In Figure [5l we consider in detail how 
the velocity in the wind direction is related to gas den- 
sity. The gas in these contours is between 2 kpc and 
50 kpc above the disk. The lower density gas is stripped 
early to a higher velocity, as shown in the top right, while 
higher density gas has almost no velocity parallel to the 
wind because the dense clouds are harder to strip and 
accelerate. This systematic correlation between density 
and velocity (differential stripping) has been observed in 
NGC 4438 with multiwavelength observations (Vollmer 



et al. 2009). There is a similar effect in the no cooling 
case (not shown), but it has a somewhat different origin: 
it is caused by the exponentially declining density dis- 
tribution in the disk. The gas density decreases towards 
the edge of the disk, which is the only part of the disk 
stripped by the wind. 

Later, as shown in the right panel of Figure O even 
higher density gas is eventually accelerated by the wind, 
although it does not reach the velocities of the lower den- 
sity gas. This gradual acceleration of denser material 
with time is not seen in the no cooling run, again because 
only the lower density gas at the edges is ever stripped. 
The differential velocity as a function of gas density is 
seen throughout the tail in all of our runs. It is also no- 
table that negative velocities, denoting gas falling back 
onto the disk, are more often seen in lower density gas 
in all of our runs. In the no cooling run, this is because 
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Fig. 5. — Contour plots of gas mass at different velocities parallel to the wind (or perpendicular to the disk) for the Tmin = 8000 K case. 
The contours are plotted as a function of gas density. Again, only gas cells with at least 25% of their mass originating from the galaxy are 
counted, and we only include gas between 2 and 50 kpc above the disk. We can see differential stripping and acceleration. See section [3?3] 
for discussion. 



there is only low density gas in the tail. In the cooling 
cases, this may be because only the lower density gas can 
be moved into the shadow of the disk by turbulence. This 
result differs from Kapferer et al. (2009), who find that 
dense H I gas falls back onto the disk. Kapferer et al. 
(2009) are not the only ones to predict gas re-accretion 
(e.g. Vollmer et al. 2001; Roediger & Briiggen 2007), 
although it has not been clearly observed. 

We can now compare our simulations to observations. 
The stripped gas moving more slowly than the ICM wind, 
at 900 km/s, is marginally consistent with (although 
somewhat larger than) the observations of Oosterloo & 
van Gorkom (2005), who find the gas from the galaxy is 
only accelerated by 550 km/s, not the 1500 km/s the au- 
thors estimate as the galaxy velocity. However, if we con- 
sider the gas velocity 120 kpc above the disk, it is closer 
to 700 km/s. Also, recall that the contours in Figured 
include gas at all densities. From Figure [5l we see that 
higher density gas moves more slowly than lower density 
gas. In fact, at any height above the disk, gas with den- 
sities greater than about 0.1 cm~^ moves more slowly 
than 900 km/s, and is frequently centered at about 500 
km/s. This puts our results into good agreement with 
the observations of Oosterloo & van Gorkom (2005). 

4. COMPARISON TO OBSERVATIONS 
4.1. HI 

In this section we consider H I column density, com- 
paring projections of our simulations to observations. We 
use recent deep observations (Chung et al. 2007; Ooster- 
loo & van Gorkom 2005) to select the minimum observ- 
able column density of H I to be 10^^ cm~^. As discussed 
in Section 12. 4[ we use Cloudy to determine the neutral 
fraction given a temperature and density and then apply 
that value in our projection routine. At the resolution 
of our simulation, we find gas with high column densi- 
ties at very large distances from the galaxy, and for a 
long time after stripping begins (see Figure [6]). The in- 
dividual H I clouds frequently have a head-tail structure, 
indicating Kelvin-Helmholtz stripping, and indeed the 
gas clumps are largely shredded in the tail, losing their 
coherency. This differs from Kapferer et al. (2009), who 
found that Kelvin-Helmholtz stripping does not affect 



the cool clouds in their stripped tails (note that Agertz 
et al. (2007) have recently argued that SPH has difficulty 
following the Kelvin-Helmholtz and Ray leight- Taylor in- 
stabilities). There is a gap between the galaxy and region 
with the most H I gas at later times. Oosterloo & van 
Gorkom (2005) also observe a gap in H I column density 
in the long tail from NGC 4388. 

Figure E] is the H I map from the T^in = 8000 K 
run. The difference between the Tmin = 8000 K and 
^min = 300 K runs (not shown) is minor, similar to the 
density distribution (see Figure [1]). The Tmin = 300 K 
run generally has a few more clouds, with slightly higher 
H I densities, where the total gas surface density is high- 
est. The gap seen in the T^in = 8000 K run is at the same 
height above the disk, but smaller in the T^in = 300 K 
run. At late times (600 Myr after the wind has hit the 
galaxy), the width of the H I tails is about the same 
in the two runs, less than 40 kpc wide, and there is a 
paucity of high surface density clouds in both runs. 

The observations in Figure [6] are at unrealistically high 
resolution, even for H I interferometric maps. We can 
also "observe" our simulation by smoothing to lower res- 
olutions, like those in typical observations; this is done 
in Figure H The first panel has a resolution of 1.2 kpc, 
similar to the observations of Chung et al. 2007. The 
middle panel has a resolution of 1.4 kpc x 7.4 kpc (to 
compare with Oosterloo & van Gorkom 2005), and the 
right panel has a resolution of 43 kpc (to compare with 
Vollmer & Huchtmeier 2007). As our smoothing size in- 
creases we find fewer clouds because some of our high 
density clouds are small enough that they are smoothed 
below our threshold column density of 10^^ cm~^. These 
three projections are all at 500 Myr after the wind has 
hit the galaxy, and are all from the Tmin = 8000 K run. 
As in Figure [6l the lowest contour is at the minimum 
column density of current observations, 10^^ cm~^. 

We begin by comparing our left panel in Figure [71 to 
the results of Chung et al. (2007), finding that our tails 
are generally longer. However, our tail is from a galaxy 
that has been stripped for 500 Myr, while the authors 
argue that the tails they observe are in the early stages 
of stripping. Indeed, about 50 Myr after the wind hits 
(see the left panel of Figure [6]) our tail is well within 
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Fig. 6. — The column density of H I gas in the Tmin = 8000 K case shown at the 38 pc resolution of our simulation 100 Myr (left), 250 
Myr (center), and 500 Myr (right) after the wind has hi t th e galaxy. The greyscale shows the log of the column density in units of cm~^. 
The longest tail is again cut off at 212 kpc. See Section ITTI for discussion. 



the range observed in these outer Virgo cluster galax- 
ies. Also, by examining their Figure 2, we see that the 
estimated ram pressures of all but one of the observed 
galaxies are less than our simulated ram pressure, so the 
tails will lengthen more slowly. Finally, our measured 
tail length is the actual length, while the tails observed 
by Chung et al. (2007) are likely to be projected at an 
angle to the line of sight. 

Next, we compare our results to the observations of 
Oosterloo & van Gorkom (2005), who observe an H I 
tail 110 X 25 kpc in size, and estimate the tail to be a 
few hundred Myr old. See the central panel in Figure 
[3 In either of our cooling cases, we can observe tails 
longer than this one 500 Myr after the wind has hit the 
galaxy. Interestingly, in the central panel, there is a gap 
with no emission, just as is seen in the tail of NGC 4388. 
The range of surface densities in the tail in this panel 
are similar to those found by Oosterloo & van Gorkom 
(2005), although our maximum H I surface density is 
about a factor of two less than they observe. In the 
^min = 300 K case at this resolution, there is a clump 
with a surface density above 3 x 10^^ cm~^, the star- 
forming lower limit, 200 kpc from the disk 500 Myr after 
the wind has hit. 

This tail of H I clouds would not have been found by 
Vollmer & Hutchmeier (2007), as we show in the right 
panel of Figure [3 At this very low resolution of about 43 
kpc, only the more distant gas in the tail is observable, 
and is more distant from the galaxy than the ~60 kpc 
search radius adopted in that paper. 



We find that the ability of the gas to radiatively cool is 
vital to the existence H I tails. As we discuss in Section 
13. 2[ as long as radiative cooling in the stripped gas is 
more efficient than adiabatic and turbulent heating, H I 
clouds will exist in our simulated tails. If heat conduc- 
tion is also not efficient enough to counteract radiative 
cooling, long H I tails consisting of dense clouds should 
be fairly common. 

In both of the cooling runs, the tail has clouds that 
are above the column density of disk gas that is associ- 
ated with star formation, 3 x 10^^ cm~^ (Schaye 2004). 
If these clouds could form stars, ram pressure stripped 
gas could result in the trails of "fireballs" observed by 
Yoshida et al. (2008) in RB 199. 

4.2. Ha 

Next, we turn to Ha emission. The minimum observ- 
able surface brightness that we use for Ha emission is 2 
X 10~^^ erg s~^ cm~^ arcsec"^ (see Sun et al. (2007) 
and references therein). Similar to our calculation of hy- 
drogen neutral fraction, we use Cloudy to determine the 
Ha emissivity given the gas temperature and density. 

We find, as shown in Figure [HI highly structured, long 
tails of Ha emission. Note that we do not include UV 
radiation from star formation or AGN (except for the 
metagalactic background, as described in section 12. 4p . 
The Ha emission roughly follows the H I column density, 
which is likely where the gas is densest and where com- 
petition between compressive and turbulent heating and 
radiative cooling is taking place, as discussed in Section 
13. 2[ The Tmin = 8000 K case shows significantly more Ha 
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Fig. 7. — The column density of H I gas observed at a resolution of 1.2 kpc (to compare with Chung et al. 2007), 1.4 kpc x 7.4 kpc (to 
compare with Oosterloo &; van Gorkom 2005), and 43 kpc (to compare with Vollmer & Hutchmeier 2007). These projections are 500 Myr 
after the wind has hit the galaxy for the standard case with cooling to 8000 K. The greyscale shows the log of the column density in units 
of cm~^. We add a red contour to the left panel, denoting a column density of 3 x lO^*-* cm~^, a lower limit for the column density of gas 
associated with star formation (Schaye 2004). 



emission than the Tmin = 300 K run, which is due to two 
reasons: first, the gas is kept close to the temperature 
where collisional excitation can result in emission. The 
other difference is that in the Tmin = 300 K simulation, 
much of the gas is at higher densities, where it is opti- 
cally thick to the metagalactic ionizing background and 
so produces little Ha emission (and it is not hot enough 
to emit due to collisional excitation). 

Turning to observations, it is important to note that 
many of the observations of Ha emission are closely 
linked to star formation or AGN activity in a galaxy 
(which we do not include in the simulations). However, 
in the case of NGC 4438 (Kenney et al. 2008), the emis- 
sion is not close to either of these ionizing sources. The 
brightness of the filaments matches well with the Ha 
emission in our simulation with radiative cooling to 8000 
K, with the lowest observed surface brightness emission 
at the level of our third contour. The overall filamentary 
structure and large gaps between emission regions in our 
projected Ha emission also matches the observations of 
Kenney et al. (2008). 

Gavazzi et al. (2001) see emission slightly above the 
brightness of our third contour, although it is difficult 



to discern the structure of the observed tails. They con- 
sider it likely that the gas was ionized within the galaxy 
and has remained ionized in the tail, although they do 
not rule out heating within the tail. We believe that tur- 
bulent and adiabatic heating of stripped cool gas causes 
much of the emission we observe in our projections, be- 
cause of the close alignment of the H I and Ha emitting 
gas added to our earlier discussion of differential acceler- 
ation (Section 13. 3p . 

As discussed in section 12. 4[ we computed our H I and 
Ha maps using Cloudy, assuming a typical cloud thick- 
ness of 100 pc, which is appropriate for our resolution. 
To explore our sensitivity to this assumption, we also 
generated maps based on cloud size of 10 pc (the 'thin' 
case), and 1000 pc (the 'thick' case). 

The thick case is very similar to our standard 100 pc 
cloud size, with slightly more H I in both of our cool- 
ing runs. In the T^in = 8000 K run, the Ha emission 
is nearly identical, possibly because collisional excitation 
dominates the emission. In the Tmin = 300 K run, the 
thicker clouds emit a little less Ha at our lowest observ- 
able brightness - although the second contour level is 
nearly identical. 
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Fig. 8. — The Ha surface brightness contours that could be 
observed using current observational depths (Sun et al. 2007). 
Again, these are for 100 Myr, 250 Myr, and 500 Myr after the 
wind has hit the galaxy. The left column is the simulation with 
radiative cooling to 8,000 K and the right side has cooling to 300 K. 
Note that the Ho; emission is filam entar y in structure and follows 
the H I column density. See sect ion |4?2] for discussion. 



The thin cloud approximation does result in lower H I 
column densities, but in both of the cooling cases there 
are still dense clouds at star forming column densities 
500 Myr after the wind hits the disk. The Ha emission 
is considerably brighter, by about a factor of 8 in the 
Tmin = 8000 K run. However, in the T^in = 300 K run, 
the emission is much brighter, frequently by a factor of 



about 16. 

Another concern is that, by using Cloudy in this way 
to generate our emissivities, we are assuming ionization 
equilibrium. Given the short cooling time, it is possi- 
ble to be out of equilibrium which could systematically 
bias the emission. To check how robust this approxima- 
tion is, we used the simulation which followed the non- 
equilibrium ionization fractions of H and He. We com- 
puted the emissivity two ways: first, with our standard 
Cloudy procedure which uses only the total density and 
temperature of the gas, and second, based on the actual 
ionization fractions predicted in the simulations based 
on a set of analytic fits to the Ha emissivity due to re- 
combination (Spitzer 1978), and the Ha emission due to 
collisional excitation (Osterbrock 1989). We found that 
both methods produced very similar Ha maps (although 
note that this run did not include cooling below T = 10^ 
K). 

4.3. X-ray 

Finally, we turn to X-ray emission. As noted in the 
introduction, only a few X-ray tails are known. We show 
predicted X-ray images in the 0.5-2 keV band in Figure[9l 
None of the contours are bright enough to be observed 
assuming an observational limit of 5 x 10~^ erg s~^ cm~^ 
(based on the luminosity and surface area of the tail in 
Sun et al. 2006). A bow shock is visible early in the 
simulations, but most of the emission comes from gas in 
the wake. Lowering the radiative cooling floor to 300 K 
does not effect the X-ray surface brightness substantially. 

There are a number of possible difference between our 
simulation and the X-ray observations of tails. First, 
since we do not include energetic feedback from super- 
novae, the gas in the tail that originated from the galaxy 
could not have been heated by super novae explosions, as 
suggested by Sun & Vikhlinin (2005). 

However, in the case of ESQ 137-001, Sun et al. (2006) 
conclude that the X-ray emission they observe is from 
mixing of the cool galactic gas with the hot ICM. This 
type of heating leading to X-ray emission could in fact 
be happening in our simulation, but just does not lead 
to observable emission. The ICM in our simulation is 
about 60% of the temperature and a third of the density 
of the ICM in A3627. Our galaxy velocity is 1.4 times 
larger than the estimated velocity of ESO 137-001, so our 
ram pressure is 63% as large as that experienced by the 
observed galaxy. However, at a distance of 200 kpc from 
the cluster center, a galaxy velocity of 2000 km/s is much 
more likely (Tonnesen & Bryan 2008), which results in 
our ram pressure being only 16% of that experienced by 
ESO 137-001. Based on their fits, the authors predict 
an X-ray tail temperature of 9x10^ K and a density of 
-|^Q-26 g cm~^. We have no tail gas from our galaxy in 
any of our simulations at that temperature and density. 
However, if we shift our contour plots in Figure [2] by the 
ratio of the ICM pressure of A3627 and our simulated 
cluster, a small amount of gas is at the required density 
and temperature. This suggests that the ICM conditions 
may be important for producing X-ray bright tails. 

5. DISCUSSION 

As we have demonstrated in the previous section, our 
stripping simulations with cooling result in very different 
wakes than in runs without cooling. This impacts both 
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Fig. 9.— The predicted X-ray surface brightness in the 0.5-2 keV band for the Tmin = 8000 K run at 100 Myr, 250 Myr and 500 Myr 
after the wind has hit. Note that none of the gas in the tail is bright enough to be observed at current observational limits. 



the morphology of the tail, the density and temperature 
structure of the gas, and the observational diagnostics. 
In this section, we begin by discussing one point which 
we have noted throughout: the differences between the 
Tmin = 8000 K and T^in = 300 K simulations. Then we 
investigate the effect of resolution and finally turn to a 
discussion of the physical mechanisms which we do not 
include in this simulation. 

5.1. Radiative Cooling Floor 

As noted earlier, we used two different values for our 
radiative cooling floor (8000 K and 300 K) in order to 
explore in a simple way the potential impact of processes 
which we do not include in the simulation. We found that 
most of our predictions are insensitive to this parameter, 
which gives us confidence in our results. However, there 
were some differences between the runs. This included 
the detailed density-temperature structure of gas in the 
tail (at the low-T end). The biggest difference was the 
amount of Ha emission, and the survival of H I clouds, 
with Tmin = 300 K predicting longer lived clouds and less 
Ha emission. 

It is unclear which Tmin value is preferable. Because 
of the lower Jeans length in the lower temperature floor 
runs, we resolve the resulting fragmentation less well in 
the run with a minimum temperature of 300 K. Although 
the Tmin = 300 K run is desirable because it allows for 
the existence of the low temperatures that are found in 
the ISM, it is not clear if this run is physically more 
realistic as we do not include some effects such as small- 
scale turbulence, cosmic rays and magnetic fields, all of 
which may provide a source of effective pressure in low 



temperature regions. 

5.2. Resolution 

We now discuss the impact of resolution on our re- 
sults. To do this, we can compare our standard simula- 
tions with runs with lower mass resolution, resulting in 
less refinement in the tail. With either refinement crite- 
ria, the gas in the disk is mostly refined to a resolution 
of 38 pc when the winds hits the disk, but in the case 
of more mass refinement, the disk fragmentation is bet- 
ter resolved, resulting in a larger number of lower mass 
clumps. This affects the evolution of the disk. In our 
radiatively cooled cases, the gas fragments from the edge 
inward. With lower mass resolution, the resulting clouds 
are more massive and can perturb the inner disk, caus- 
ing fragmentation and the formation of over- and under- 
dense regions. In contrast, in the highly resolved cases 
that we discuss in this paper, the disk has fewer holes in 
the central disk, and therefore stripping does not occur 
in the central 7 kpc. Also, stripping occurs more slowly, 
although a slightly larger amount of gas is stripped from 
the galaxy (this is very similar to our discussion regard- 
ing the case with cooling to 300 K in TB09). 

With the lower mass resolution run, most of the gas in 
the tail is refined to either 304 pc or 152 pc, and none of 
the tail is refined to 38 pc. This is in contrast to the more 
refined tail, in which most of the gas is refined to either 
152 pc or 76 pc, and there are a number of small areas 
refined to 38 pc. We go through our results discussed in 
this paper, highlighting the similarities and differences of 
different resolution runs. 

The morphology of the stripped gas differs early in 
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the stripping process, but by 500 Myr after the wind 
has hit the galaxy, the stripped gas looks very similar. 
In Figure [10] we show the most different projection of 
either of our comparison sets of runs (Tmin = 8000 K 
and Tmin = 300 K), using the T^in = 8000 K cases at 250 
Myr after the wind has hit the galaxy. The lower mass 
resolution case has fewer dense clouds at the edges of the 
disk, so most of that gas is quickly stripped. However, 
the higher resolution run has a slightly longer tail because 
the very low density gas in that run has been accelerated 
more quickly to larger distances. The galaxies with more 
refinement have more and smaller clouds in the tail, as 
expected. 

The gas density and temperature distributions are very 
similar. Although the gas distribution remains similar, 
at late times the runs with lower resolution in the tail 
have less high density gas and a slightly lower maximum 
density. This is because clouds do not condense to as 
small sizes in the lower resolution tail. 

The tail velocity profiles are also very similar, with sim- 
ilar widths and the same maximum velocities reached at 
the same distance from the galactic disk. All radiatively 
cooled runs have a consistent velocity width over time. 
There is also some fallback onto the disk in the lower res- 
olution run, because larger gas clouds in the disk cause 
fragmentation in the central regions and therefore also 
result in gas stripping in those regions. This centrally 
stripped gas can then be protected by disk gas that is 
rotating or that moves to fill the center of the disk (dis- 
cussed in detail in TB09), and will be shielded from the 
stripping wind. This protection allows even more dense 
gas to have slightly negative velocities in the lower resolu- 
tion cases. Thus, the gas punching through the bottom of 
the disk in our lower resolution runs and in the Kapferer 
et al. (2009) simulations may be resolution dependent. 

The H I maps are somewhat different, again because of 
the increased fragmentation in our highly resolved tails. 
Compare the two projections in Figure [11] to the first 
panel in Figure [T] We find that the Tmin = 300 K run 
produces more H I in our maps than the T^in = 8000 K 
run, a more marked difference than in our standard high 
mass resolution runs. This is because when we use the 
more stringent mass criterium to resolve the tail, only 
the lower cooling fioor creates high density clouds that 
are above our refinement threshold and therefore more 
highly resolved. The H I tails are generally shorter, and 
the Tmin = 8000 K run has much less gas at large dis- 
tances in the tail. Also, both cooling fioors have tails 
that are more connected to the disk in the lower resolu- 
tion runs. This occurs for two reasons: first, the tail is 
shorter because the gas is less dense and is more easily 
heated by the turbulence in the tail. Second, gas can be 
stripped from the central region of the disk and then pro- 
tected by rotating or infiowing disk gas. Therefore, gas 
that is immediately behind the disk may be protected 
from the wind and can survive there for longer. 

Ra emission is largely unaffected, because the same 
type of turbulence exists in the different resolution cases 
despite the difference in size of the dense clouds. The 
interaction of the slower moving dense clouds with the 
faster low density wind still causes small scale turbu- 
lence that can heat the nearby gas. The difference in 
resolution does not change the brightness of the emis- 
sion because the temperatures are the same, resulting in 




Fig. 10. — Contour plots of total gas surface density, comparing 
two Tmin = 8000 K cases with different mass resolution 250 Myr 
after the wind has hit. The left figure is our standard high resolu- 
tion run, and the right is the lower mass resolution run. Because 
there is less gas clumping in the disk, especially in the outer re- 
gions, the lower resolution run has a smaller remaining disk. There 
is more fallback in the lower resolution run, resulting in some of the 
gas punching back through the disk. Although the dense clouds in 
either case are close to the same distance above the disk, the less 
dense gas has clearly been accelerated more quickly from the more 
highly refined disk. 



similar amounts of collisional excitation. If the difference 
in clouds sizes were enough to greatly affect the optical 
depth of the clouds. Ha emission would decrease. The 
main difference is the distribution of emission, which fol- 
lows the H I clouds in the tail in any run. Therefore, as 
with the H I gas, the Ha emission is closer to the disk. As 
we discussed above, the cooling floor has a much larger 
effect. 

The X-ray emission is affected very little. As we dis- 
cussed in Section 14. 3[ the surrounding ICM pressure 
likely needs to increase for X-ray emission to be observ- 
able. 

In summary, the different results caused by using dif- 
ferent resolutions can be attributed to two main causes. 
First, and most importantly, less resolution in the tail 
results in less fragmentation to smaller dense clouds (al- 
though there is more fragmentation in the T^in = 300 
K run than in the T^in = 8000 K run). Second, gas 
stripped from the galactic center can be protected by 
shadows from the wind by the remaining disk gas, and 
fall back towards the disk. 

5.3. Heat Conduction 

Heat conduction could be important for the survival 
of cool clouds in the ICM and for Ha emission. If heat 
conduction is an efficient way to transport heat from the 
ICM to cold, stripped gas, then the survival time of H I 
clouds would be less than predicted in this paper, and the 
length of the tails would be shorter. The fact that the 
length of our tails matches observations of either shorter 
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Fig. 11. — Contour plots of H I gas column density at 1.2 kpc 
resolution, comparing the two lower mass resolution runs with ra- 
diative cooling. These are both projected 500 Myr after the wind 
has hit the galaxy, and can be compared to the left panel of Fig- 
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: 8000 K is on the left, and 



300 K is on the 



right. Notice that there is more dense gas closer to the disk. As 
we discuss in Section 15.21 centrally stripped gas is able to sit in 
the protected lee of the disk, and dense gas does not survive as 
far from the disk because it does not collapse into as small dense 
clouds. There is more dense gas in the Tmin = 300 K case because 
the lower cooling floor results in smaller clouds. 



tails by Chung et al. (2007) for our T^^ 



8000 K 



case, or long tails like that observed by Oosterloo & van 
Gorkom (2005) using the Tmin = 300 K simulation, ar- 
gues against a strong role for heat conduction, at least 
for the current conditions. 

We can perform an estimate of the evaporation time 
for a typical cloud if heat conduction is not suppressed. 
We follow Cowie & McKee (1977), as in Vollmer et al. 
(2001). In our clouds, we find that the mean free path 
for ions is comparable to or greater than the tempera- 



ture scale length, so we need to use the saturated heat 
flux equations. Solving for the evaporation time, we find 
that using a cloud radius of 100 pc, and a number den- 
sity of about 0.1 cm , t^yap 6 X 10 years. However, 
this time can be lengthened if there are magnetic fields 
perpendicular to the temperature gradient (Cowie & Mc- 
Kee 1977; Cox 1979). Vollmer et al (2001) argue that 
magnetic fields could increase the evaporation time by 
nearly an order of magnitude, although the actual value 
is highly uncertain. Although McKee & Cowie (1977) 
predict that a cloud in our modeled ICM would be evap- 
orated before being radiatively cooled, the cooling time 
for a cloud of density 0.1 cm~^ and temperature 10^ K 
is 10^ years. A higher density shortens the cooling time 
and lengthens the evaporation time, so radiative cooling 
may also slow down the effects of heat conduction. 

5.4. Magnetic Fields 

Our simulation does not include magnetic fields or cos- 
mic rays. In general in the Milky Way, the magnetic field 
is measured to be about 2-3 /iC (Men et al 2008). How- 
ever, there are measurements of the magnetic field in 
molecular clouds (Crutcher 1991 and references therein) 
of up to 10^ /iC. If some of the denser clouds we see in 
the tails in our simulations had strong magnetic fields, 
these fields could protect the clouds from ablation while 
they move through the wind, in addition to lowering the 
efficiency of heat conduction. We do not know to what 
extent magnetic fields could counteract turbulent heat- 
ing and ablation of stripped clouds because we do not 
know the amount of structural support magnetic fields 
would give stripped clouds. 

The first long stripped tails were observed in radio con- 
tinuum, with stretched magnetic field lines downstream 
of the galaxy (Gavazzi et al. 1995). Including magnetic 
fields in our simulation would likely also result in radio 
continuum emission (assuming a cosmic ray population 
was also stripped), although because our ram pressure is 
less than the ram pressure estimated to be affecting the 
three galaxies in A 1367, the magnetic fields would be 
stretched less. 

5.5. Star Formation 

A molecular cloud may survive between a few and a 
few tens of Myr (Blitz & Shu 1980; Larson 2003; Hart- 
mann 2003). Because we do not include star formation, 
our dense clouds of gas do not evolve into less dense gas 
that could be stripped. If our dense regions of gas could 
turn into low density regions through star formation {p < 
10-23 g QYR-^^^ Q^j. galaxies that include radiative cool- 
ing could be entirely stripped of gas over the lifetime of 
the dense clouds. This could result in a tail without a 
gap between the majority of the stripped gas and the 
galactic disk. 

We find that overdensities in the stripped tails can 
survive to large distances, however, our projections also 
show that they are being ablated by the ICM wind. Our 
next step is to add star formation, which could affect 
stripping in a number of ways. Ram pressure could in- 
duce star formation in the disk, using up disk gas and 
resulting in lower mass tails and fewer dense clouds in 
the tail. If there are dense clouds in the tail, we will dis- 
cover whether stars can form in stripped gas, or if Kelvin- 
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Helmholtz stripping will destroy clouds before they can 
collapse. 

6. CONCLUSIONS 

We have run detailed galaxy simulations including ra- 
diative cooling to understand the morphology of and 
emission from gas tails stripped by ram pressure in a 
cluster environment. We compare three cases: a run 
without cooling, one with cooling to a minimum temper- 
ature floor of Tmin = 8000 K to approximately account 
for non-thermal pressure support, and one to Tmin = 300 
K. Our main conclusions are as follows: 

1. Including radiative cooling creates a clumpy ISM, 
which results in faster stripping of lower density gas, and 
the production of longer tails. A clumpy disk also has 
more holes for the ICM wind to pass through, which re- 
sults in significantly less fiaring in the stripped material, 
and narrower tails in the cooling runs, in better agree- 
ment with observed tails (e.g. Chung et al. 2007). 

2. We find that the clumpy ISM also produces a much 
wider range of densities and temperatures in the wake, 
resulting in a very different morphology for the tail. Al- 
though the stripped gas is mostly at the ICM pressure, 
it shows a wider distribution in pressure, as well as a 
systematic deviation from the constant pressure line for 
temperatures below T = 10^ K, due to radiative cooling 
operating more quickly than compressive heating. 

3. The stripped gas in the cases with radiative cooling 
does not accelerate to the ICM wind speed. This is likely 
because the dense clumps are harder to accelerate than a 
smooth distribution of gas. Also, the gas is differentially 
stripped and differentially accelerated through the whole 
length of the tail for the duration of the simulation. 

4. We compute realistic H I maps of our tails and 
compare to H I observations of stripped tails. With either 
a radiative cooling fioor of 8000 K or one of 300 K we find 
H I tails with long lifetimes (more than 500 Myr), and 
large lengths (well over 100 kpc). When we consider our 
run with cooling to 300 K, we find that the tails are longer 
and clouds can survive at high density for a longer time. 
Most tails that have been observed in Virgo have shorter 
lengths, like our tails at very early times, however, many 
of these tails seem to be in the early stages of stripping 
(as in Chung et al. 2007), and may have only experienced 
a short amount of time at large ram pressure (see the 
orbits in Vollmer et al. 2001). The few longer tails that 
have been observed are well matched by our results at 
late times in the cooling runs. 

5. We can also map the Ha emission from our runs, 
and find filamentary structures much like those found in 
deep observations of stripped tails (e.g. Sun et al. 2007). 
Both the Tmin = 8000 K and the T^in = 300 K case 
emit in Ha, but the simulation with the higher cooling 
fioor has more emission and a higher maximum surface 
brightness in the tail. The filamentary structure in either 
radiative cooling run reproduces the observations of the 
Ra tail of NGC 4438 (Kenney et al. 2008), while the 
^min = 8000 K run has some emission that is bright 
enough to be observed at the Kenney et al. (2008) limit. 
This Ha emitting gas is near the H I gas, indicating 
that it is created in the tail rather than surviving as Ha 
emitting gas stripped from the disk. 

6. None of the simulated tails have X-ray surface 
brightnesses high enough to match observations of X-ray 



tails (note that some observed H I or Ha tails also do not 
show X-ray emission) . This may either be because we do 
not preheat the galactic gas with super novae explosions, 
or that our ICM is not hot and dense enough for mixing 
to cause bright X-ray emission. 

It is clear from our results that adding radiative cool- 
ing greatly improves the realism of simulated tails. This 
is seen by comparing both the morphology and velocity 
structure of the simulated tails to observations. We find, 
as in observations, that there are only two types of bright 
emission from our tail: in our case H I and Ha emission. 
It remains to be discovered whether all three types of 
emission can be at observable levels in a single tail or 
if the environment necessary for X-ray emission would 
destroy or heat the cooler clouds. However, we can use 
our results and a few analytic estimates to make an ini- 
tial guess: if our ICM had the same pressure as that of 
A3627 we would likely have some gas in our tail with the 
necessary density and temperature for detectable X-ray 
emission (see Section 14. 3p . At the same time, it may 
be that heating from compression or turbulence would 
still not be enough to heat clouds faster than they would 
be radiatively cooled, so we might still see H I in the 
tail. Finally, Ha emission should be present as long as 
there is some intermediate density gas (near T ~ 10^ 
K), so it may not be greatly affected by the ICM pres- 
sure. To complete our comparison to the case of ESO 
137-001, we have a lower ram pressure than that affect- 
ing the observed galaxy, and Sun et al. (2007) find that 
the ram pressure is likely to have triggered star forma- 
tion throughout the disk. This may have heated the cool 
gas clouds surrounding the star forming regions so much 
less cool gas was left to be stripped. As we discuss in 
Section 15. 5[ we intend to add a star formation prescrip- 
tion to a future simulation, which will allow for a better 
understanding of how star formation can affect the gas 
tail. 

The H I column density projections indicate that 
higher resolution observations (similar to that of Chung 
et al. 2007) of cluster galaxies would find more massive 
tails. The fact that the H I tails can have large gaps 
between the higher column density regions means that 
there could be tails that would not have been found by 
Vollmer & Hutchmeier (2007), and points to a need for 
large surveys at the highest possible resolution. 

Our simulations predict that only slightly deeper obser- 
vations of tails in Ha should reap great rewards of much 
more structure. It is interesting that although less dense 
gas moves more quickly than the dense gas throughout 
the tail, the emission from Ha is nearly perfectly aligned 
with the H I column density, indicating that gas in our 
simulations is heated to Ha emitting temperatures in- 
situ, rather than being stripped while ionized and sur- 
viving for hundreds of Myr. 

We can consider what observational signatures of ram 
pressure stripping our simulations predict are likely. We 
find dense H I clouds can survive far from the disk and 
for a significant amount of time. This means that in 
order to determine that any H I structure is not sur- 
viving stripped gas, other information such as kinematic 
data indicating rotation, is necessary (as in Haynes et 
al. 2007). In addition. Ha emission can be observably 
bright even without in-situ star formation. However, the 
high column densities of some of the H I clouds leads us 
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to believe that stripped gas can condense to form stars 
that could be formed far from the galaxy, indicating that 
even a stellar component in a stripped tail is not enough 
to prove that a gravitational interaction took place. 

The high resolution combined with the addition of ra- 
diative cooling in our simulation allows us to follow the 
interaction of dense cool clouds in the tail with the hot 
gas of the ICM. Unlike some previous work, we find ab- 
lation and destruction of the clouds are important, as 
they move within the wind, limiting their survival to less 
than found by Kapferer et al. (2009). It is important 
to reiterate that because it is difficult to make a physi- 
cal case for either cooling fioor, we present both 8000 K 
and 300 K cooling fioor s. Although many of the main 
results about the tail attributes are very similar between 



the two cases, they do result in meaningfully different 
observable profiles in Ha. Finally, although these simu- 
lations are arguably the most realistic stripping simula- 
tions performed to date, we have discussed in detail the 
limitations involved in our work in ^ 
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